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Executive  Summary 

Characteristics  of  ground  clutter  in  stationary  and  airborne  radar  systems  is  studied  herein. 
Both  the  monostatic  and  bistatic  configuration  are  considered.  The  focus  is  on  how  the 
clutter  characteristics  are  affected  by  the  manner  in  which  observation  is  performed.  The 
factors  considered  include  the  transmitted  waveform,  pulse  compression,  radar  platform 
motion,  antenna  pattern,  and  system  geometry. 

The  analysis  starts  with  a  brief  review  of  the  space-time  processing.  This  review  also 
shows  the  manner  in  which  clutter  estimation  is  involved  in  the  space-time  processing.  The 
basic  model  for  the  clutter  analysis  is  then  presented.  It  is  followed  by  the  derivation  of  the 
clutter  autocorrelation  function  for  three  cases,  namely,  clutter  observed  from  a  stationary 
monostatic  radar,  from  an  airborne  monostatic  radar,  and  from  an  airborne  bistatic  radar. 

It  should  be  pointed  out  that,  in  this  study,  no  consideration  is  given  to  the  variation 
in  clutter  strength  due  to  the  different  incidence  and  observation  angles.  However,  this 
problem  cannot  be  solved  until  a  constant  Gamma  model  is  developed,  and  verified  through 
experiment,  for  the  out-of-plane  clutter  in  bistatic  configuration.  The  assumption  about  the 
scatterer  spectrum  may  also  be  relaxed  to  allow  for  the  derivation  of  clutter  process  from 
basic  principles,  especially  for  the  derivation  of  non-Gaussian  clutter. 

The  goal  of  this  effort  is  two-fold.  First,  through  the  analysis  of  the  effects  of  the  system 
parameters  on  the  clutter  characteristics,  criteria  for  secondary  data  selection  and  for  signal 
design  in  space-time  processing  can  be  developed.  Another  application  of  this  analysis  is  that 
it  provides  as  a  basis  for  the  development  of  the  space-time  processing  for  bistatic/multistatic 
radars.  A  preliminary  study  is  presented  in  this  report  to  show  how  the  multistatic  case 
departs  from  its  monostatic  counterpart.  The  procedure  for  further  development  is  also 
discussed. 
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Chapter  1 
Introduction 


Detection  of  slowly  moving  objects  close  to  the  ground  or  sea  surface  by  means  of  an  airborne 
radar  is  seriously  impaired  by  clutter  returns.  These  returns  exhibit  Doppler  shift  which 
varies  with  the  velocity  of  the  radar  platform  and  the  direction  of  the  antenna  beam.  In 
addition,  even  absolutely  stationary  clutter  has  a  Doppler  spectrum  due  to  the  fact  that 
reflections  of  scatterers  from  different  parts  of  the  beam  have  different  Doppler  shifts. 

Clutter  cancellation  is  effectively  achieved  by  filtering  out  the  portion  of  the  Doppler 
spectrum  occupied  by  the  clutter  return.  However,  the  target  may  also  be  eliminated  in 
the  process  if  Doppler  is  the  only  condition  used  to  distinguish  the  clutter  and  the  target, 
since  the  sidelobe  clutter  may  overlap  the  target.  This  problem  is  encountered  by  the  tradi¬ 
tional  Moving  Target  Indicator  (MTI).  In  contrast,  adaptive  space-time  processing  further 
distinguishes  the  clutter  and  the  target  return  in  the  angle  dimension  by  employing  arrays 
of  sensors.  Thus,  by  assigning  appropriate  weights  to  the  sensors,  the  notch  of  the  clutter 
cancellation  filter  is  placed  along  the  line  where  the  clutter  lies  in  the  angle- velocity  plane. 

To  select  the  weights  to  the  sensors,  it  was  shown  that  the  weights  which  maximize  the 
signal-to-clutter-plus-noise  ratio  (SCNR)  are  related  to  the  clutter  covariance.  In  practice, 
however,  the  clutter  covariance  is  substituted  with  its  estimate  when  forming  the  weights, 
since  this  quantity  is  generally  not  known  a  priori.  This  estimate  is  obtained  from  clutter 
samples  in  real-time.  As  a  result,  space-time  processing  greatly  improves  the  output  signal- 
to- noise  ratio  with  minimal  a  priori  information  about  the  changing  environment. 

Numerous  research  efforts  from  both  the  radar  and  the  sonar  communities  have  been 
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devoted  to  the  analysis  of  space-time  processing.  Development  of  this  scheme  still  lacks  in 
two  areas,  however.  First,  there  are  few  studies  for  the  bistatic/multistatic  sj stems  in  the 
open  literature.  This  type  of  system  configuration  consists  of  one  or  multiple  transmitters 
and  receivers  that  are  separated  in  space.  The  advantage  is  the  decreased  vulnerability 
over  its  monostatic  counterpart;  nevertheless,  the  increased  complexity  in  system  geometry 
introduces  difficulties  for  effective  realization  and  analysis. 

Another  area  of  the  space-time  processing  scheme  that  requires  further  investigation  is 
how  the  performance  is  affected  by  the  clutter  statistical  estimation.  This  question  is  of 
concern  especially  when  spatial  averaging  is  involved  in  the  clutter  estimation,  since  clutter 
from  different  locations  may  not  have  the  same  spectrum. 

In  this  research  effort  we  consider  the  aforementioned  two  problems,  but  the  emphasis  is 
on  the  study  of  the  problems  associated  with  clutter  estimation.  To  set  the  stage  for  later 
development,  we  give  in  the  following  a  brief  review  of  how  the  space-time  processing  was 
derived  for  the  monostatic  system.  In  Section  2  we  describe  the  assumptions  and  method¬ 
ology  used  in  analyzing  the  clutter  statistical  characteristics.  In  Chapter  2,  we  study  the 
clutter  characteristics  as  observed  from  a  stationary  monostatic  system,  an  airborne  mono- 
static  system,  and  an  airborne  bistatic  system.  In  Chapter  3,  we  examine  how  the  Doppler 
spread  in  the  ground  clutter  affects  the  clutter  spectrum  and  clutter  correlation  function. 
A  preliminary  study  of  space-time  adaptive  processing  for  bistatic/multistatic  is  given  in 
Chapter  4.  Chapter  5  discusses  the  results. 

1.1  Space-time  processing  in  monostatic  radar  system 

In  this  section  we  briefly  review  the  manner  in  which  space-time  processing  was  derived  for 
the  monostatic  radar  system.  The  treatment  follows  closely  to  that  given  in  [1]. 

Filtering  of  the  received  signal  to  improve  the  signal-to-clutter-plus-noise  ratio  is  essential 
to  the  detection  of  a  moving  target  close  to  the  ground  or  sea.  However,  schemes  which 
involve  processing  only  in  the  time  or  the  frequency  domain  are  ineffective;  they  may  remove 
both  the  clutter  and  the  target  when  the  sidelobe  clutter  exhibits  the  same  Doppler  as  the 
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target.  This  situation  is  illustrated  in  Figure  1.1,  where  v  and  V  are  the  speed  of  the  radar 
platform  and  of  the  target,  respectively,  and  6  are  the  azimuth  and  the  depression  angles, 
respectively,  and  A  is  the  wavelength.  For  simplicity,  we  assume  the  depression  angle  is  small. 

In  most  cases,  the  clutter  and  the  target  can  be  further  distinguished  in  the  angular 
dimension,  as  depicted  in  Figure  1.2.  This  figure  shows  that  the  added  degree  of  freedom 
in  clutter  processing  enables  one  to  avoid  the  problem  of  rejecting  both  the  clutter  and  the 
target.  Clearly,  the  angular  information  can  be  obtained  when  an  array  of  sensors  is  used. 

Suppose  that  the  echo  from  a  target  is  received  by  a  set  of  N  sensors.  Let  each  sensor 
yield  I\  temporal  samples  of  the  received  signal  such  that  a  total  of  L  =  N  x  A'  samples  are 
obtained  at  the  receiver.  Further  assume  that  the  sampling  is  performed  at  the  I/Q  channel 
output.  Denote  the  entire  sample  data  set  by  the  column  vector  s  =  [si  •  •  •  sl]t,  where 
Sfc.A-  =  1.2. ••  •  ,L  is  the  kth  signal  sample  and  the  superscript  T  denotes  the  transpose 
operation.  Similarly,  denote  the  noise  vector  by  nc  =  [nj  n2  •  •  •  Further,  let  M„e 

denote  the  noise  covariance  matrix.  The  detection  of  the  presence  of  a  target  is  by  putting 
the  received  vector  through  a  filter  with  weights  w  =  [u>i  u*2  •  •  •  wi]T. 

The  output  signal-to-noise  ratio  of  the  filter  is  given  by 

wHs\ 

d2  =  — hTtw  —  (1*1) 

wHMncw 

where  H  denotes  the  Hermitian  operation.  Using  the  Schwarz  inequality,  it  can  be  shown 
[1]  that  <P  is  maximized  by 

w  =  kMl)s  (1.2) 

where  Jt  is  a  constant  not  equal  to  zero.  It  can  also  be  shown  [1]  that  the  weights  given  by 
(1.2)  maximize  the  detection  probability,  at  a  given  false  alarm  rate,  in  a  Gaussian  noise 
environment. 

Realization  of  (1.2)  requires  the  knowledge  of  M„c.  Since  M„e  is  seldom  known  a  priori, 
it  needs  to  be  estimated.  Denote  the  estimate  of  M„c  by  M„e.  Thus,  the  weights  to  the 
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Figure  1.1:  Clutter  Doppler  spectrum  for  airborne  monostatic  radar. 
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Figure  1.2:  Clutter  in  the  angle-Doppler  space 
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sensors  are  actually  calculated  from 


w 


(1.3) 


There  are  two  methods  commonly  used  for  the  estimation  of  the  clutter  covariance  matrix. 
One  is  to  use  the  sample  covariance  matrix,  which  is  the  maximum  likelihood  estimate  for 
the  Gaussian  clutter.  Another  approach  is  the  maximum  entropy  method  (MEM),  which 
involves  the  estimation  of  the  clutter  power  spectrum.  With  either  method,  spatial  averaging 
across  several  range  cells  is  generally  to  increase  the  available  samples,  under  the  assumption 
that  the  clutter  is  homogeneous  in  a  close  neighborhood.  However,  clutter  samples  obtained 
from  different  range  cells  correspond  to  clutter  return  from  different  range,  while  clutter 
returns  from  different  range  may  not  have  the  same  Doppler  spectrum  due  to  the  factors 
previously  described.  The  geometry  effects  on  the  matched  filter  output  due  to  a  point  target 
was  analyzed  in  [2].  In  the  next  chapter  we  will  consider  the  factors  that  affect  the  clutter 
statistical  characteristics. 

Another  part  of  the  realization  of  (1.2)  is  the  choice  of  s.  This  choice  has  the  effect  of 
canceling  out  the  undesired  interference  and  is  called  a  steering  signal.  To  illustrate  the  use 
of  steering  signals,  consider  the  following.  Assume  the  array  is  linear  with  sensor  element 
spacing  /.  Then,  the  relative  phase  of  the  kth  sensor  due  to  a  signal  arriving  at  an  angle  tp 
with  respect  to  the  array  normal  is  given  by 


Vk  = 


sinyj,  k  =  1,2,---,.V 


(1.4) 


The  above  relationship  indicates  that  signal  arriving  in  different  directions  induces  different 
relative  phases  among  the  sensors  .  This  feature  is  exploited  by  adjusting  the  phases  of  s 
to  tune  out  the  undesired  interference.  In  Chapter  5  we  present  a  preliminary  study  of  how 
the  above  analysis  can  be  extended  to  the  multistatic  case. 


1.2  Basic  model  for  clutter  analysis 

In  this  section  we  discuss  the  assumptions  and  methodology  commonly  adopted  in  the  lit¬ 
erature  on  the  analysis  of  clutter.  This  consideration  provides  us  with  a  basis  for  the  devel- 
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opment  of  a  statistical  description  of  the  clutter,  to  be  used  in  our  study  of  the  estimation 
performance  and  the  effects  on  the  space-time  processing. 

The  most  adopted  assumption  about  clutter  is  that  the  clutter  is  a  collection  of  indepen¬ 
dent  point  scatterers.  The  number  and  reflectivity  of  the  scatterers  are  then  modeled  as 
random  numbers.  The  appropriate  assignment  of  the  underlying  probability  distribution  for 
different  types  of  terrain  have  been  the  subject  of  many  research  efforts.  The  occurrence  of 
a  scatterer  in  time  is  generally  assumed  to  be  governed  by  Poisson  process,  e.g.,  in  [3].  This 
assumption  leads  to  Gaussian  clutter  by  further  assuming  that  the  number  of  scatterers  is 
large.  Other  models  for  the  number  of  the  scatterers  were  also  proposed  in  the  literature. 
In  particular,  a  recent  paper  [4]  shows  that  clutter  of  Spherically  Invariant  family  can  be 
obtained  by  using  Poisson  mixture  distribution. 

There  are  system  factors  which  affect  clutter  characteristics.  Clearly,  the  ground  area 
illuminated  depends  on  the  transmitting  and  receiving  apertures  and  the  system  geometry. 
The  size  of  a  resolution  cell  also  determines  how  many  scatterers  contribute  to  the  clutter 
return  [5].  Furthermore,  the  clutter  spectrum  is  a  function  of  the  relative  motion  between 
the  scatterers  and  the  radar  platform.  Finally,  there  is  the  effect  of  the  relative  position  of 
the  scatterers  and  the  radar  platform,  which  affects  the  strength  of  the  reflections,  as  well  as 
the  location  and  shape  of  the  clutter  spectrum.  A  rather  thorough  treatment  on  the  effects 
of  these  factors  can  be  found  in  Middleton’s  work  [6]- [10]. 

Our  analysis  will  differ  from  the  previous  work  in  two  aspects:  we  will  include  pulse 
compression  and  give  an  explicit  description  of  the  geometry  relationships.  In  modern  radars, 
pulse  compression  is  used  to  improve  the  resolution  capability.  Pulse  compression  is  achieved 
through  raising  the  bandwidth  of  the  signal  by  modulation,  at  transmission,  and  matched 
filtering,  at  reception.  As  a  result,  the  samples  available,  i.e.,  the  output  of  each  range 
resolution  cell,  are  the  result  of  matched  filtering.  Thus,  it  is  of  interest  to  know  how 
the  estimated  clutter  characteristics  are  subjected  to  change  due  to  matched  filtering,  as 
compared  to  using  samples  from  sampling  the  received  signal  directly. 

Another  issue  which  is  generally  overlooked  is  an  explicit  derivation  of  the  geometry  re¬ 
lationships.  Rather  than  giving  a  set  of  equations  to  describe  the  geometry,  we  will  solve 
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them  explicitly  and  discuss  the  implications.  This  will  lend  some  insight  as  how  system 
performance  varies  with  geometry. 

Specifically,  we  will  assume  that  the  transmitted  signal  is  a  periodic  pulse  sequence.  Re¬ 
flection  from  the  ground  is  modeled  as  a  doubly-spread  return  [11].  This  is  a  modification 
to  KlemnTs  approach  used  in  [12],  where  the  clutter  from  one  range  cell  is  modeled  as  a 
Doppler-spread  return.  One  recent  work  adopted  a  similar  approach  [13].  The  clutter  re¬ 
turn  is  assumed  to  be  the  superposition  of  reflections  from  a  number  of  independent  point 
scatterers.  At  any  instant,  the  region  from  where  the  scatterers  contribute  to  the  clatter 
return  corresponds  to  the  extension  of  one  range  resolution  cell,  which  is  obtained  from  pulse 
compression.  The  spatial  distribution  of  the  scatterers  is  governed  by  a  probability  density 
function.  The  number,  n,  and  the  reflectivity,  g,  of  the  scatterers  are  modeled  as  random 
numbers.  The  strength  of  the  reflections  also  depend  on  the  antenna  transmit  pattern  Gj 
and  receive  pattern  Gr,  the  transmitted  waveform  s(f ),  and  the  distance  the  signal  traveled. 

We  further  have  the  following  assumptions.  First,  the  terrain  is  assumed  to  be  flat  [14]. 
The  range  of  depression  angle  of  interest  is  assumed  small  such  that  the  range  resolution 
is  determined  by  the  signal  pulse  duration  [15].  Range  ambiguities  are  not  considered  at 
present  analysis.  The  probability  density  function  describing  the  location  of  the  scatterers 
does  not  change  within  the  observation  interval.  The  fluctuation  of  the  reflectivity  of  the 
scatterers  during  the  observation  interval  can  be  modeled  as  a  wide-sense  stationary  process. 

The  analysis  will  concentrate  on  how  the  clutter  characteristics  vary  across  range  cells,  i.e., 
a  quantified  description  of  the  change  across  range  cells.  We  will  identify  how  the  observed 
clutter  characteristics  are  subject  to  change  due  to  the  transmitted  waveform,  the  manner  in 
which  the  received  signal  is  processed,  platform  motion,  different  shapes  of  antenna  pattern, 
and  system  geometry.  Finally,  the  analysis  will  cover  both  the  monostatic  and  bistatic  radar 
systems. 
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Chapter  2 

Observation  of  Clutter  Characteristics 


In  this  chapter  we  will  study  how  the  clutter  characteristics  vary  as  a  function  of  the  system 
parameters,  including  the  transmitted  signal,  pulse  compression,  radar  platform  motion, 
antenna  pattern,  and  system  geometry.  The  analysis  starts  with  an  examination  of  the  clutter 
characteristics  when  observed  from  a  stationary  monostatic  radar  platform,  and  arrives  at  the 
generalized  case  where  the  radar  platform  is  airborne  and  bistatic.  The  derivation  is  based  on 
the  point-scattering  model,  by  which  the  clutter  return  is  written  as  a  linear  superposition 
of  the  reflections  from  point  scatterers.  The  secondary  scattering  effect  is  assumed  to  be 
negligible:  the  occurrence  of  the  scatterers  is  then  assumed  to  be  governed  by  the  spatial 
Poisson  process. 


2.1  Clutter  observed  from  a  stationary  monostatic  sys¬ 
tem 

In  the  following  we  analyze  the  clutter  characteristics  when  it  is  observed  from  a  monostatic 
radar  which  is  monostatic  and  is  mounted  on  top  of  a  stationary  tower.  The  derivation  is 
based  on  the  point-scattering  model  where  the  clutter  return  is  written  as  a  linear  superpo¬ 
sition  of  the  reflections  from  point  scatterers  distributed  on  the  ground. 

In  this  system,  position  of  a  point  on  the  ground  can  be  specified,  relative  to  the  radar 
site,  by  giving  the  height  k  of  the  platform  and  the  azimuth  angle  <p  and  distance  R  of  the 


9 


(2.1) 


point.  The  projection  of  R  on  the  ground  is  called  the  ground  distance  Rg, 

Rg  =  VR2  -  h2. 

For  convenience,  we  assume  that  the  distance  is  measured  in  terms  of  the  wave  propagation 
time.  The  depression  angle  6 ,  measured  from  the  local  horizontal  of  the  platform,  positive 
downward,  can  be  obtained  from 

6  =  cos"1  %.  (2.2) 

R 

The  scenario  considered  and  the  relationships  of  the  parameters  are  illustrated  in  Figure  2.1. 

The  clutter  return  is  assumed  to  be  a  linear  superposition  of  the  reflections  from  point 
scatterers  distributed  throughout  the  illuminated  area.  Denote  by  Dr  the  region  swept 
by  the  transmitted  pulse  from  the  minimum  to  the  maximum  detection  range.  Using  the 
subscript  k  to  denote  the  quantities  attributed  to  the  kth  scatterer,  the  total  reflection  from 

the  scatterers  in  Dr  can  be  written  as 

Nd  z 

rc{t)  =  ^^( ^)Gr(v?i:,  Ok)s{t  -  2 Rk),  (2.3) 

*=o  Kk 

where  A 'd  is  the  total  number  of  scatterers  in  Dr.  Let  td  be  the  signal  round-trip  travel  time 
to  some  point  in  the  range  resolution  cell  Dc  under  consideration:  also,  denote  the  correlation 
operation  by  0.  The  clutter  return  observed  in  Dc  is  given  by 

nc(td)  =  ^GT{(Pk,0k)GR((pk,0k)s(^  ~  0  s(t  —  td).  (2.4) 

Denoting  the  auto-correlation  function  of  s(t)  by  />*(•), 

Ps(t)  =  /  s(t)s*(t  - r)dt,  (2.5) 

J —oo 

the  expression  (2.4)  can  be  written  as 

nc(td)  =  5Z  0k)Gfi(<pk,  9k)p,(td  -  2Rk)-  (2.6) 

Jfc=0  Kk 

The  clutter  correlation  is  defined  as 

Mnc(tut2)  =  E[nc(ti)n;(t2)].  (2.7) 
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The  evaluation  of  (2.7)  requires  the  joint  probability  distribution  of  the  number,  the  location, 
and  the  reflectivity  of  the  scatterers.  For  this  purpose,  we  have  the  following  assumptions. 
First,  the  occurrence  of  the  scatterers  is  assumed  to  be  governed  by  a  spatial  Poisson  process 
[16], [17].  Due  to  the  fact  that  an  illuminated  area  usually  consists  of  a  number  of  regions 
having  different  scattering  nature,  e.g.,  patches  of  sand,  vegetation,  and  artificial  objects, 
we  assume  that  the  illuminated  area  can  be  divided  into  Nr  such  regions  and  the  intensity 
function  of  the  Poisson  process,  fic(Rg,ip),  can  be  written  as 

t*c{Rg, *p)  =  Me,.'  for  {Rg,  <p)  €  A„  *  =  1,  •  •  • ,  Nr,  (2.8) 

where  .4,  is  the  ith  region  in  Dr.  It  follows  that  the  scatterers  in  each  region  are  uniformly 
distributed  when  conditioned  on  a  fixed  number  of  scatterers  [16],  with  the  mean  value 
//CI|A,j.  where  |A,j  is  the  size  of  A,-.  Further,  we  assume  that  the  scatterers  are  statistically 
independent  of  each  other,  the  number  of  the  scatterers  is  independent  of  the  scatterer 
reflectivity  and,  in  each  A„  there  is  a  correlation  function  />5.,(*)  describing  the  fluctuation 
of  the  reflectivity  of  the  scatterers,  i.e., 

E  [gk(ti)gi  (<2)]  =  PgA1  2  “  *1  )hi,  (2-9) 

for  the  A-th  and  /th  scatterers  in  A,-. 

There  is  one  more  issue  to  be  addressed  before  we  can  substitute  (2.6)  into  (2.7)  to  evaluate 
the  clutter  correlation.  The  problem  arises  when  the  time  axis  of  n(ti)  and  n(f2)  has  different 
origins.  For  example,  this  is  the  case  when  n(tx)  and  n{t2)  are  obtained  from  different  pulses 
of  a  transmitted  pulse  train.  Let  the  difference  in  time  origins  be  Td,  with  n(t2)  having  a 
later  time  origin.  Then,  denoting  by  D *  and  Dr9  the  domains  of  the  depression  angle  and 
the  ground  distance  of  the  scatterers,  respectively,  and  by  Nc  the  number  of  A,-  that  fall 
within  Dc.  we  have 

=  !tw,Ah  -  I.)  L  t  |Cr(v’,e)CW*>.«)|! 

I=1  JURg,i 

-jjiP.it  2  -Td-  2  R)p:(U  -  2  R)dRgdg>.  (2.10) 
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In  the  above  expression,  the  integration  with  respect  to  Ra  can  be  changed  to  with  respect 
to  R  by  using  (2.2),  i.e., 

dRg=sec0dR.  (2.11) 

Also,  the  integration  interval  Dc  can  be  extended  to  infinity,  since  the  product  ps(t2  -Td  — 
2R)p"s(ti  -  2 R)  is  negligible  outside  of  Dc ;  then, 

^  ,  ,00 

-lfnc(tl,  ^2)  =  r  PcSPQ.iidi  ~~  ^l)  /  l  Q)Gr{y'Q)\2 

»=1  JDv..i  J 00 

sec# 

■-RA-Psih  -I'd-  2  R)P;{t1  -  2  R)dRd<p.  (2.12) 


Assume  that  the  approximation  of  the  factor  secO/R*  with  a  constant  sec  Sd/M,, 


Od  —  cos' 


Rd 


(2.13) 


within  Dc  does  not  result  in  appreciable  error.  Then, 

sec  Op  Na 


Mnc(tut2)=  \  j" 

nD  ,=l  J<x> 

■p,{t2  -Td-  2 R)pl(U  -  2 R)dRdip.  (2.14) 


Also,  assume  that  the  antenna  patterns  are  approximately  constant  with  respect  to  R  within 
Dc.  such  that 

\Gt{<p,0)\  =  |Gr(r,,^z?)|  (2.15) 

and 

I  Gr(<p,  0)|  =  |Gfl(^,  Od)  I  •  (2.16) 

It  follows  that 

sec  On  (-d  t 

Mnc(h,t2)  =  p4  ]C PcjPg.ifa  ~ii)  /  \Gt(v>,Od)Gr((p,0d)\2 

kd  ,=1  Jdv,< 

•  l°°  Ps{i2  -Td-  2R)P;(h  -  2 R)dRd<p.  (2.17) 

J  OO 

Since  the  integral  involving  the  two  auto- correlation  functions  in  (2.17)  can  be  written  as 
r  Ps{t2  -Td-  2 R)p](t1  -  2 R)dR  =  2  {p,{t2  -  h  -  Td)  ®  p,(t2  -  tx  -  Td)} ,  (2.18) 

J  -'OO 
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we  have 


sec  Or)  r- t 

Mnc(ti,ii)  =  0D4  'HvcjPgAh-t i)  /  |Gr(^> Od)Gr{ip, Od)\2 dip 

ltrD  i=i  •/£\» 

•  W<2  -  <1  -  Td)p's(t2  -  <1  -  Trf)}  .  (2.19) 

The  clutter  power  spectrum  ,Mnc(-)  is  the  Fourier  transform  of  A/,,,..  Denote  the  Fourier 
transforms  of  p5(-)  and  s(-)  by  Ps(-)  and  S(-).  respectively.  From  the  fact  that 

|2 


and 


we  have 


/CO  I  y-OO 

{/9,(x)  0  p4(x)}  e~j2*yxdt  =  I  ps(x)e~j2iryxdt 

-GO  |  J  —  OO 

r  p,(x)e-^dx  =  \S(y)f , 

J — OO 


(2.20) 


(2.21) 


Mnc (/)  =  ^E^fa.*/)  ®  L  \Gt{^d)Gr{^0d)U<p.  (2.22) 

2.2  Clutter  observed  from  an  airborne  monostatic  sys¬ 
tem 


In  this  section  we  generalize  the  previous  analysis  by  assuming  the  radar  platform  is  airborne. 
Through  this  examination,  the  effects  of  the  relative  motion  between  the  radar  platform  and 
the  illuminated  scatterers  on  the  observed  clutter  characteristics  will  be  clear. 

The  relative  motion  between  the  radar  platform  and  the  scatterers  induces  Doppler  shift 
v.  Explicitly,  v  is  a  function  of  the  platform  speed  v,  the  direction  of  v,  the  transmitted 
wavelength  A,  and  the  location  of  the  scatterer.  Let  (<pv,0v)  be  the  direction  of  v;  then, 

2v 

v  —  —  sin(v?  -  (pv)  cos (0  -  0V).  (2.23) 

Comparing  with  (2.3),  where  the  radar  platform  was  assumed  to  be  stationary,  the  reflections 
from  the  scatterers  in  Dr  becomes 

rc(0  =  E  TkGT((pk,0k)GR(ipk>0k)s(t  —  2 Rk)ej27”/k(t~2R^ dRg.  (2.24) 

*= 0  Kk 


14 


Suppose  the  return  from  range  ring  Dc  is  obtained  by  correlating  the  received  signal  with 
s(t  -  #rf)exp[/27r/u(<  -  <</)],  where  td  is  the  travel  time  to  some  point  in  Dc  and  /D  is  the 
Doppler  shift  with  respect  to  a  point  in  Dc.  Consequently,  the  observed  clutter  return  is 
given  by 

nc(td)  =  |  ^  -  2JRfc)<r»'2,r^^-2,l*>  J  ®  «(*  - 

Denote  the  time-frequency  correlation  function  of  s(t)  by  <1>4(-,  •), 

*.(r.  /)  =  ^  (t  +  0  s*  (t  -  0  e~*z}tdL  (2.26) 

The  expression  (2.25)  becomes 

Nd 

nc{U)  =  Y.  %GT(<Pk,0k)GR(vk,0k)<t>s(td  -  2 Rk,fD  -  (2.27) 

k= 0  Kk 

With  the  same  assumptions  and  procedures  used  in  deriving  (2.17),  we  have 

sec  0  n  t 

Mnc(t  1,<2)  =  —^4 — Yl  Hc'ipg.iit*  ~  tl)  I  \Gt(<P,0d)Gr(<P,0d)\2 

i=i  JDf,i 

•  f°°  $4(*2  -  7*  -  2R,  fD  -  i/)$;(tj  -  2/2,  /D  -  v) 

J  OO 

.tM!D-^h-h-Ta)dRd(p  (2.28) 

Due  to  the  dependence  of  v  on  R,  the  right-hand  side  of  (2.28)  is  not  a  function  of  (t2  —  h) 
only:  this  indicates  a  non-stat ionary  process.  Then,  it  is  of  interest  to  find  out  if  there  are 
cases  where  the  nonstationary  effect  is  not  significant. 

Consider  the  effect  of  a  small  deviation  A/?  of  Rq  on  0.  The  corresponding  change  of  0, 
A0,  can  be  found  from  the  relationship 


h  =  (R  +  A  R)  sin(0  +  A0)  =  R  sin  0. 
Using  the  small  angle  approximation 

sin  A0  =  A0 


and  straightforward  manipulation  of  (2.29),  we  have 

-A  R/R 


A0  =  tan  0 


1  +  AR/R‘ 


(2.29) 


(2.30) 


(2.31) 
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In  practice,  it  is  usually  true  that  the  width  of  the  range  ring  is  much  smaller  than  the 
minimum  detection  range;  in  other  words,  A R/R  <C  1,  for  all  R.  Thus, 


A  0  =  tan  9 


-A  R 
R  ' 


With  the  approximation 


cos  AO  =  1, 

the  value  of  the  cosine  factor  around  (Od  —  0V)  can  be  written  as 


(2.32) 

(2.33) 


cos{0d  +  A 9-  9V)  =  cos (Od  -  0V)  —  s\ti(0d  -  6V)A0.  (2.34) 


From  (2.23),  (2.32)  and  (2.34),  the  Doppler  shift  in  the  reflections  from  the  scatterers 
located  near  Rd  can  be  approximated  by 


v  =  ^  sin(<^>  —  <pv)  cos (Od  —  9V)  +  ^  sin(</?  —  <pv)  sin(0p  —  0V)  tan  Od-tJ--  (2.35) 
A  A  Kd 


It  can  be  seen  from  (2.35)  that,  by  denoting  the  first  term  on  the  right-hand  side  by  vp  and 
the  second  term  by  Aj/;  the  quantity  Av  can  be  ignored  when 


sin(#£)  —  0V)  tan  0qAR/ Rd 
cos  (Od  —  Ov) 


=  tan  (Od  —  0V)  tan  @D 


A  R 
Rd 


<  1. 


(2.36) 


Therefore,  it  is  sufficient  that 

tan(0£>  —  0V )  tan#£>  <  1  (2.37) 


for  (2.36)  to  be  true,  e.g.,  Od  <  tt/4  is  sufficient  when  0V  =  0.  The  clutter  correlation  (2.28) 
becomes 


t2) 


sec  0  d 
Rd 


Na 

^c,i ~  tl) 
i=l 


•  f  \GT(^jD)GR(^M\2^MfD'UD){t2~ti~Td) 

JDv,i 

too 

•  /  9.(t2-Td-2R,fD-v)$;(U-2RyfD-v) 

Joo 


(2.38) 


In  the  rest  of  this  analysis,  we  consider  only  the  case  where  Av  <C  vd- 
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In  a  manner  similar  to  the  derivation  of  (2.19),  we  have 

nD  i= l 

•  /  |Gr(^^)Gft(v?,^)|2ei!r(/D-,/D)(t2~‘,-:r<') 

•  {$s(t2  -h-  Ti ,  /o  -  t/)$;(f2  -  <1  -  2*  fD  -  v)}  d<p.  (2.39) 

The  clutter  power  spectrum  is  readily  obtained  by  taking  the  Fourier  transform  of  the  right- 
hand  side  of  (2.39).  From  the  fact  that 

/“  {*.(!,»)  9  *.(!•.»)}  e-»’w‘ix  =  ||_”  .  (2.40) 

and 

/“  *.(i,  y)ritmix  =  5  (t»  +  |)  5-  (w  -  |)  ,  (2.41) 

we  have 

AM/)  =  £  Pc*  L  I Gt(v,  0d)Gr(<p,  0d) I2 

*H>  i= i 

■  In,'  (/  -  ®  |S(/)S  (/  -  (/d  -  "d))|2J  if.  (2.42) 

2.3  Clutter  observed  from  an  airborne  bistatic  system 

In  the  following  we  consider  the  case  where  the  radar  platform  is  airborne  and  bistatic.  The 
effects  of  geometry  on  clutter  will  be  examined  in  detail. 

The  Cartesian  coordinate  system  is  used  to  describe  the  three  dimensional  bistatic  geome¬ 
try.  Let  the  transmitter  site  be  the  origin;  this  choice  is  arbitrary  and  will  be  discussed  later. 
Let  the  local  horizontal  plane  of  the  transmitter  be  the  x-y  plane.  The  vertical  line  passing 
through  the  transmitter  is  chosen  to  be  the  2-axis,  and  the  intersection  of  the  x-y  plane 
and  the  vertical  plane  containing  the  transmitter  and  the  receiver  is  chosen  to  be  the  x-axis; 
y-axis  is  defined  according  to  the  right-hand  rule.  Further,  assume  that  the  depression  angle 
is  measured  from  the  x-y  plane,  positive  downward,  and  the  azimuth  angle  is  measured  from 
y-axis,  positive  clockwise. 
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For  simplicity  of  expression,  measurements  relative  both  to  the  transmitter  and  to  the 
receiver  are  used.  In  particular,  we  assume  that  the  angle  measurements  relative  to  the 
receiver  are  obtained  by  translating  the  origin  from  the  the  transmitter  site  to  the  receiver 
site.  The  same  symbols  for  the  angle  and  distance  measurements  used  in  the  monostatic 
case  are  kept  for  the  present  analysis;  however,  subscripts  T  and  R  are  added  to  indicate 
that  the  quantity  is  measured  relative  to  the  transmitter  or  the  receiver.  For  example,  Rj 
is  the  distance  of  a  point  from  the  transmitter  and  Rjt  is  the  ground  projection  of  Rj. 

Rt9  =  \JW^T  (2.43) 


The  symbol  R  is  now  used  to  denote  the  total  distance  the  signal  traveled:  R  =  Rj  +  Rr. 
Finally,  let  L  denote  the  distance  between  the  transmitter  and  the  receiver.  The  scenario 
considered,  together  with  the  relationships  of  the  parameters,  is  illustrated  in  Figure  2-2. 
Specifically,  the  quantities  Rr ,  <pr  and  Or  can  be  expressed  as  functions  of  Rt,  and  Oj 
by  knowing  the  relative  positions  between  the  transmitter  and  the  receiver.  Let  depression 
angles  of  the  receiver  with  respective  to  the  transmitter  be  Ort',  clearly,  the  azimuth  angle 
of  the  receiver  relative  to  the  transmitter  is  ?r/2.  From  the  law  of  cosines, 

Rr  =  ]^3  +  L 2  cos2  Ort  -  2Rt<iL  cos  0RT  sin  y> t  +  •  (2-44) 


Using  the  relationships 


Rt  sin  Oj  =  Rr  sin  Or  +  L  sin  Ort 


and 


Rt  cos  <?t  cos  Qt  =  Rr  cos  cos  Or, 


it  can  be  shown  that 

Or  =  sin-1  ((i?r  sin  Qt  —  L  sin  Ort)/Rr] 

and 


<Pr  =  cos-1  {( Rt  cos  y>r  cos  Qt  \R?r  cos2  Qt  +  L2  cos2  Ort 
•  +2 RtL  cos  Qt  cos  Ort  cos(y5 t  —  v’rt)]-1^2}  • 


(2.45) 


(2.46) 


(2.47) 


(2.48) 
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Figure  2.2:  Geometry  relationships  for  the  airborne  bistatic  radar 
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In  comparison  with  the  clutter  model  (2.25),  the  factor  R~ 2  representing  the  propagation 
loss  becomes  (RjRr)~1-  Also,  the  Doppler  shift  depends  on  the  bistatic  geometry.  Let  the 
transmitter's  speed  relative  to  the  ground  be  vj,  in  the  direction  (<pvr,0v T),  and  the  receiver’s 
speed  relative  to  the  ground  be  vr,  in  the  direction  (<pVR,0vR)-  The  resulting  Doppler  shift 
is  given  by 

l*T  l?p 

v  =  —  sin(v?r  -  VvT)  cos{0T  -  0VT)  +  -j-  sin(y>/j  -  <pVR)  cos{0R  -  0VR).  (2.49) 

Hence,  the  clutter  return  is  written  as 

MO  =  |e  Z?  %--GT(vTk-0Tk)GR(<pRk,9Rk)s(t  -  /?fc)ej2,r‘'*^_H*)| 

®s(t  -  tD)ei2*Mt-iD)  (2.50) 

Then,  similar  to  the  monostatic  case,  we  assume  that  the  reflectivity  of  the  scatterers  are 
statistically  independent  from  the  number  and  from  the  location  of  the  scatterers,  and  the 
occurrence  of  the  scatterers  is  governed  by  the  non-homogeneous  spatial  Poisson  process. 
The  bistatic  clutter  correlation  function  becomes 

•Y*  r  roo  1 

MnAU-h)  =  J^f*c,iPg,i(t2-tl)  /  \Gt{PT,9t)Gr((Pr,0r)\  —  —  - 

1=1  JD*T.i  J~‘x’  ( RtRr ) 

■$s(t2  -Td-R,fD-  v)$ms(ti  -  RJd-  v) 

-t>-T<>)dRTad<pT.  (2.51) 

It  should  be  noted  that,  in  (2.51),  the  integration  over  Rt9  is  due  to  the  choice  of  the 
transmitter  site  as  the  origin  of  the  coordinate  system;  the  integration  over  Rjg  can  be 
changed  to  over  R.  Define  the  factors 

A  =  L  cos  0rt  sin  p>j  (2.52) 


and 

B  =  hj-  —  h2R  —  L2  cos2  0rt.  (2.53) 

Using  the  fact  that  Rr  =  R—  Rt  and  the  relationship  (2.44),  it  can  be  shown  that 

Rt,  =  [A(#2  +  B)  +  Ry/{R*  +  Bf  -  4 K^{R2  -  A2)]/  [2 (R2  -  A2)] ,  (2.54) 
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for  R  >  '  r  +  \J L2  cos2  6rt  +  and 


fiTg  =  [a(#2  +  B)  ±  ^(^FO2-^^2-^2)]/  [2(i?2  _  A2^  (2.55) 

We  have  the  following  observation  for  the  case  R  <  hj  +  \J L2  cos2  0rT  +  The  region 

R  <  hr  -  y[£ cos2  6rt  +  h\  on  the  ground  is  inside  the  contour  defined  by  the  intersection 
of  the  ellipsoid  R  =  h?  +  yj L2  cos2  Or t  +  h\  and  the  ground.  When  the  region  R  <  hj  + 
\J L1  cos-  9rt  +  h2R  is  larger  than  R  >  h.R  +  \J L2  cos2  6rt  +  hj-,  it  is  more  convenient  to  set 
the  origiz  of  the  coordinate  system  at  the  receiver  instead  of  the  transmitter.  In  such  regions, 
there  are  two  intersections  with  the  constant  range  contour  in  the  transmitter’s  line  of  sight 
and  may  result  in  ambiguities;  we  only  consider  the  case  where  R  >  hT  +  \JL2  cos2  0rt  -f  h2R. 
The  analysis  can  easily  be  extended  to  the  case  where  R  <  hj+y/L2  cos2  9rt  +  h2R  by  adding 
the  additional  clutter  return.  It  can  be  shown  that 

^  =  2RtRr  [(fl2  +  B)7  -  4h7T(R7  -  A2)]-1/2 .  (2.36) 

More  discussions  on  Rt,  Rr  and  Rt9  are  given  in  Appendix  A. 

From  the  foregoing  derivation,  the  clutter  correlation  function  (2.51)  can  be  written  as 

.!/„,< -  *1)  f  (  | Gt(<Pt,0t)Gr{<pr,0r)\2 

i=l  JD*T,i  J-oo 

"(jirk)2*,fe  ~Ti~R's°- ")'*'*'('>  -  a/o  - 

•2 RtRr  [(R7  +  B)7  -  4hj(R7  -  A7)]  ~,/J  dRdfr.  (2.57) 

For  simplicity  of  notation,  define  the  quantity 

Q  =  2RtRr  [(£2  +  B)2  -  4 h\{R2  -  A2)]  "1/2 .  (2.58) 

Thus, 

Mnt'-t\,h)  =  ]T)/*c,«7>a,.(*2  -  <l)  f  f  \Gt(<PT,6t)Gr(<Pr,6r)\2  \2 

i=l  JD*T-i  J~°°  \K-tKr) 

•<M*2  -Td-RJD-  »/)$;(<!  -RJd-  v)j*U°-')te-t'-TAdRdvT(  2.59) 

Suppose  Dc  is  short  enough  that  a  constant  Rr  can  be  used  in  calculating  the  factor 
Q/(RtRr )2  without  causing  appreciable  error.  Let  Rtd,Rrd  and  Qd  denote  the  quantities 
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obtained  by  substituting  Rq  into  (A.l)  and  (A. 3)  of  Appendix  A  and  (2.58),  respectively, 


e.g.. 

Qd  =  2 RtdRso  [(«£>  +  Bf  -4Af.(ff|,  (2.60) 

Then. 


Na 

*^ne(^  11^2)  =  He.ipQ.i 

t'=l 


Qp  r 

(RtdRrd)2  J-X 


|Cr(^r,  Qt)Gr(<pr,  Or)  |2 


•$a(/2  —  Td  —  R,  fD  —  v)$;(h  -  RJd  -  vW*UD-l'Ht'-t'-T^dRdipT{ 2.61) 


Similar  to  the  monostatic  case,  we  assume  that  the  antenna  patterns  are  insensitive  to  the 
change  of  R  within  Dc\  let  0td  and  0rd  denote  the  depression  angles  evaluated  at  a  constant 
Rd  in  Z)c,  e.g.. 

9Td  =  sin’1  (2.62) 

It  follows  that 


MnAU-h)  =  YtVc'iPgAt 2-U)  f  ,p  \Gt{yT,0T)Gr(<Pr,0r)\2 

t=l  *'^*7'**  \**'Tp*£Rd) 

•/  $s(h-Td-  /2,/z? -i/)$*(<i  -  R.fD-v) 

J- 00 

.tMfD-»W*-u-Td)dRdipT  (2.63) 


The  expression  (2.63)  indicates  a  nonstationary  process  since  v  varies  with  R  such  that 
(2.63)  is  not  a  function  of  (f2  —  ti).  The  condition  under  which  the  Doppler  shift  can  be 
considered  approximately  constant  with  respect  to  R  within  one  range  ring  was  analyzed  for 
the  monostatic  case;  the  derivation  can  be  extended  to  the  bistatic  case.  It  can  be  seen  from 
the  derivation  of  (2.31)  to  (2.36)  that  it  is  sufficient  to  have 

(2.64) 

and 

!«.(»«„-«„)  ton  (2.65) 

For  v  to  be  insensitive  to  the  change  of  R  in  one  range  ring  consider  the  factors  ARt/Rtd 
and  A Rr/Rrd.  It  is  expected  that  those  two  factors  are  also  much  less  than  unity  when 
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v 


A R/R  <  1  for  one  range  ring  for  fairly  general  geometry.  The  change  of  RT  due  to  a  small 
change  of  R  is  governed  by  the  differentiation  of  Rj  with  respect  to  R.  Using  chain  rule,  we 
have 


dRj  dRj  dRT 


From  (2.43)  and  (2.56),  we  have 
dRj 


dR  dRj  dR 


dR 


=  2RT9Rr  [(/?2  +  B)2  -  4 h2T(R2  -  A2)] 


-1/2 


It  follows  that 

=  2 j&Rr  [(R2  +  B)2  -  4 h2T(R2  -  A2)] _,/2  dR. 
Substituting  (A. 3)  from  Appendix  A  for  Rr ,  we  have 


dRj 

Rt 


Rt, 

Rj 


l-2A2/R2  -  B/R2 


A/R 


dR 
R  ' 


(2.66) 


(2.67) 


(2.68) 


(2.69) 


.(1  -  A2 / R2 ) ^(  1  +  B/R2)2  -  4/4//?2(l  -  A2/#2)  1  -  A2/R2 

Assume  that  the  square  of  the  ratios  hT/R,  hR/R  and  L/R  are  much  smaller  than  unity; 
it  follows  that  A2/R2  and  B/R 2  are  negligible  compared  to  unity  and  Rt9/Rt  is  close  to 
unity.  Consequently, 

dRx  /,  A\  dR 

~r7  =  v  -  r)  ~r'  ^2,70) 

From  the  definition  of  A  given  by  (2.52),  we  have  \A/R\  <  1.  Then,  we  conclude  from  (2.70) 
that  ARj/Rt  <  1  when  A  R/R  <  1  within  one  range  ring.  The  analysis  of  the  condition 
for  A RRg/ Rr  to  be  much  smaller  than  unity  follows  similar  procedure  as  that  for  ARt9/Rt: 
it  is  given  in  Appendix  A. 

The  factors  |  tan  6td  |  and  |  tan  6Rd  \  are  smaller  than  unity  by  the  previous  assumptions: 
therefore,  it  is  sufficient  to  have  |tan(0rD  -  9VT)\  and  |tan(^  -  0VT)\  to  be  smaller  than 
unity  for  (2.64)  and  (2.65)  to  be  true.  Assume  that  |tan(0j^  —  9vr)\  and  |  tan(0jD  —  6VT)\ 
are  smaller  than  unity  in  the  analysis  to  follow.  Let  i/d  denote  the  quantity 


i/d  =  -j  sin(v?r  -  w)  cos {6Td  -0vt)  +  y  sin(^fl  -  <pVR)  cos (9Rd  -  9Vr).  (2.71) 

\\  ith  the  above  analysis,  it  is  seen  that  u  is  insensitive  to  the  change  of  R  within  one  range 

ring  when  the  assumptions  axe  true,  such  that 

Na  ,  Q 

Mndh-h)  =  I>c.,7V(*2  (i?T  Br  )2  \^'r^f>T‘'^TD)GR{ifR^9RD)f 
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(2.72) 


•  {$s(*2  -  <1  -  Ti.fo  -  I'd)  ®  $m,(h  -  h  -  TdJo  -  I'd)} 

.ej”(fD~l'D)(h-ti-Td)(JxpT 


It  follows  that 


.W„(/) 


E  /•«  /D>  1  {n,  (/  -  ®  |5(/)5  (/  -  (/o  -  "o))!2} 

•  7 p  ^  |GfT(<KT,^T£,)G!fi(iJ5/i,^fiD)|2<fv?r  (2.73) 


2.4  Detection  of  target  in  clutter 

The  effects  of  ground  clutter  return  on  the  detection  performance  in  both  the  space-time 
and  time  domain  processing  are  considered  in  the  follwoing.  In  particular,  we  examine  the 
SCNR  at  the  receiver  output.  This  quantity  is  also  proportional  to  the  detection  probability 
in  the  Gaussian  interference  case. 

The  receiver  output  SCNR  was  given  by  (1.1),  and  its  maximum  is  attained  by  using  the 
weights  specified  by  (1.2).  As  a  result. 


£  =  sHR~is.  (2.74) 

Suppose  that  the  thermal  noise  has  power  o\  and  denote  the  clutter  covariance  matrix  by 
Rc:  the  expression  (2.74)  becomes 

(2.75) 


where  I  denotes  the  identity  matrix  and  Rc  denotes  the  clutter  correlation  matrix.  Specifi¬ 
cally,  Rr.  can  be  written  as 

o\l  +  -Reo  Rc, 


Rn  = 


K 


+  K 


<*/+*»  j 


(2.76) 


where  each  Rc,  is  the  correlation  matrix  due  to  the  clutter  component  of  the  output  of  the 
array  sensors  at  i  pulses  apart,  and  N  is  the  number  of  the  received  pulses.  In  other  words, 
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the  ( A\  /)th  element  of  Rc,  is  given  by  £[nCfc(f)n*  (t  —  »T)J.  It  follows  that  Rn  is  not  of  full 
rank  if  on  =  0  and  E[sk(t)sf(t  —  t'T)]  =  E[sp{t)s*(t  —  iT)].  When  this  problem  arises,  it  can 
be  avoided  by  adding  a  small  quantity  to  the  diagnal  of  Rn.  In  the  following,  without  making 
further  assumptions  regarding  the  observed  clutter  characteristics  across  the  array  sensors, 
we  will  only  consider  the  SCNR  at  the  output  of  a  single  sensor  on  a  bistatic  receiver.  A 
more  thorough  investigation  of  the  utility  of  (2.39)  and  (2.72),  especially  in  the  context  of 
space-time  processing,  is  the  focus  of  the  continuing  effort. 

The  derivation  of  the  SCNR  at  the  output  of  a  single  sensor  receiver  is  discussed  in  detail 


in  [3]: 


SCNR 


=/: 


|S«(/)I! 


rtf. 


(2.77) 


■»  <  +  M'JU) 

where  Si{j)  is  the  desired  signal  from  the  receiver  output. 

The  result  (2.77)  has  two  important  consequences.  First,  it  can  be  seen  that  the  perfor¬ 
mance  of  the  optimum  receiver  is  the  same  as  the  matched  filter  when  Ad*c(/)  is  displaced 
from  $Af)-  Therefore,  it  indicates  that  the  receiver  structure  can  be  simplified  with  properly 
designed  transmitted  signal.  This  subject  was  studied  in  detail  in  the  literature,  e.g.  [5]; 
however,  no  consideration  has  been  given  in  the  context  of  space-time  processing. 

The  second  issue  associated  with  (2.77)  is  when  S(f)  is  wide  but  £*(/)  is  narrow,  in  terms 
of  the  quantity  fo  —  vd,  i-e.,  the  Doppler  spread  in  one  resolution  cell.  If  we  can  further 
assume  that  the  illuminated  ground  clutter  is  homogeneous,  it  can  be  seen  from  (2.73)  that 
the  area  of  the  resolution  cell  becomes  the  dominant  factor  which  affects  the  receiver  output 
SC’NR.  Specifically,  when  the  the  illuminated  ground  clutter  is  homogeneous,  we  can  write 

/d  ~  vd' 


MnAf)  =  n.jDr  \v,  (f  -  0  lS(/)S(f  -  Up  -  I'd))!2! 


Qd 


(RtdRrd)2 


I  Otd)Gr(<pr,  0rd )|2  dipr 


(2.78) 


From  the  above  expression  it  is  seen  that,  by  letting  an  impulse  function  to  approximate 
Tg(f)  and  letting  a  constant  to  approximate  S(f)  in  the  correlation  calculation,  the  integral 
is  affected  by  the  geometry  factors  alone.  In  the  following  we  give  an  example  to  illustrate 
this  observation. 
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Example 


In  this  example  we  we  consider  the  clutter  return  in  a  bistatic  system  where  both  the 
transmitter  and  the  receiver  depression  angles  are  small.  We  assume  that  the  transmit 
beam  is  omnidirectional  but  the  receive  beam  is  narrow  with  a  3  dB  beamwidth  of  2  deg; 
for  simplicity,  we  further  assume  that  the  receive  beam  pattern  is  uniform  within  the  3  dB 
beamwidth  and  negligible  elsewhere.  The  above  assumption  about  the  receive  beam  is  one  of 
the  two  models  which  are  generally  adopted  in  the  literature,  while  Gaussian  beam  pattern 
is  the  other  one.  Thus,  we  have 


/p  \Gt{<Pt,  0td)Gr('Pr,9rd)\2  dtpr 

(  K'TdK'Rd  ) 


L 


Qd 


db  (RtdRrd)2 


d<pT, 


(2.79) 


where  Db  is  the  ground  intersection  of  the  beam  and  the  distance  corresponding  to  one 
range  cell.  The  quantity  on  the  right-hand  side  of  (2.79)  is  evaluated  along  three  range 
rings,  at  the  receiver  look  angle  from  -n  to  7r,  while  each  range  ring  is  specified  in  terms  of 
its  eccentricity  e,  R  =  L/t.  Figure  2.3  shows  the  geometry  of  e  =  0.4  by  having  the  baseline 
to  be  50km  and  the  range  ring,  the  inner  circle  in  the  plot,  to  be  125km.  The  results  are 
plotted  by  normalizing  with  respect  to  the  value  obtained  at  y>R  =  jt/2.  Thus,  the  plots  can 
be  used  for  any  geometry  irregardless  of  the  actual  values  of  R  and  L.  Figures  2.4  to  2.6 
give  the  results  for  e  =  0.4,  e  =  0.1  and  e  =  0.025.  respectively. 

From  these  figures  it  is  seen  that  the  geometry  factors  have  mild  to  minimal  effects  on  the 
clutter  power.  This  conclusion  is  consistent  with  results  in  the  literature  where  the  change 
of  illuminated  area  is  considered  to  be  the  dominant  factor  in  the  variation  of  the  clutter 
power.  However,  it  should  be  emphasized  that  the  applicability  of  the  results  shown  in  this 
example  relies  on  the  vailidity  of  the  simplification  that  leads  to  (2.79). 
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Figure  2.3:  Isorange  Contour  at  eccentricity  of  0.4 
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Normalized  Geometry  Factor 


Eccentricity  e=0.4 


Figure  2.4:  Normalized  geometry  factor  at  eccentricity  of  0.4 
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1 


Eccentricry  e=0.1 


Figure  2.5:  Normalized  geometry  factor  at  eccentricity  of  0.1 
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Normalized  Geometry  Factor 


Eccentricity  e=0.025 


Figure  2.6:  Normalized  geometry  factor  at  eccentricity  of  0.025 
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Chapter  3 

Doppler  Spread  of  the  Ground 
Clutter  in  Airborne  Bistatic  Radar 


This  chapter  concerns  the  Doppler  spread  of  the  ground  clutter  in  airborne  bistatic  radars. 
The  subject  of  Doppler  spread  of  the  clutter  in  the  bistatic  plane  was  treated  in  [16].  Here 
we  extend  the  development  to  three-dimension,  assuming  flat  earth.  This  assumption  is 
useful  when  the  illuminated  region  on  the  ground  is  small  such  that  it  can  be  approximated 
by  a  plane.  To  facilitate  the  analysis,  we  first  derive  some  results  on  the  three-dimensional 
bistatic  geometry.  The  ground  clutter  Doppler  spread  is  studied  in  Section  2. 

3.1  Constant  range  and  constant  Doppler  contours 

In  this  section,  we  show  how  the  intersection  of  the  constant  range  contour  and  the  ground 
can  be  described  for  an  airborne  bistatic  radar,  with  the  assumption  of  flat  earth.  The 
illuminated  portion  of  this  intersection  gives  the  cross-range  dimension  of  the  resolution  cell 
on  the  ground. 

Expressed  in  terms  of  the  Cartesian  coordinate  system,  the  constant  range  contour  in  an 
airborne  bistatic  radar  is  an  ellipsoid  given  by 

(x2  -f  y2  +  z2)  1  +  [(x  —  Lcos  0rt)2  +  y2  +  (x  -  +  ^t)2]  '  =  R-  (3.1) 

Assuming  flat  earth,  the  intersection  of  this  ellipsoid  with  the  ground  is  an  ellipse,  which  can 
be  found  by  setting  z  =  —hj  in  (3.1).  With  this  substitution  and  straightforward  algebraic 


31 


(3.2) 


manipulations,  the  expression  for  the  ellipse  is  obtained  as 

(•r-<*)2  ,  y2  . 


where  the  major  axis  is 

1  - 


R 

a=2 


(hR  -  hT)2 


b2 


ii/2  r 


R 2  — *  L2  cos2  9rt 


(hR  +  hTy 


1/2 


the  minor  axis  is 


(7?2  -  I2  cos2  ^at)1/2  L  _  (fefl  —  /ij) 


^2-I2cos2^J 

( hR  +  hT )2  11/2 


1 1/2  r 


2  [~  R2  —  L 2  cos2 

the  translation  in  the  x-axis  is 

L  cos  0R j  ( ,  h?R  —  hj 


1  - 


R2  —  L2  cos2  0Rt 


d  = 


(- 


=)• 


(3.3) 


(3.4) 


(3.5) 


2  V  R 2  —  ^2  cos2  J 

and  sqrtR2  —  L2  cos2  0Rt  >  +  /ij)  is  assumed  to  insure  that  the  intersection  exists. 


Using  the  definition  c  =  \/a2  —  62,  the  location  of  the  foci  of  the  ellipse  can  be  written  as 
( -c  +  (L  0,  —  hr)  and  (c  +  d,  0,  —  hj). 

To  determine  the  extent  of  a  range  cell,  the  coordinate  of  the  points  on  the  contour  (3.2) 
as  function  of  the  transmitter  or  receiver  azimuth  angle  is  needed.  Clearly,  the  points  on 
the  ground  at  a  given  <pj  satisfy  y  =  tan  <pxx.  Substituting  this  relationship  into  (3.2),  we 
obtain  the  coordinates 

b2d cos2  +  <*6  cos  <px\J (a2  ~  d2)  sin2  tpx  +  b2  cos2  ~pj 


x  = 


a2  sin2  (pT  +  b2  cos2  ip? 


(3.6) 


and 


b2d sin  <px  cos  PT  +  ab  tan  <p j  cos  < px\J (a2  —  d2)  sin2  ^pj  +  V2  cos2 
^  a2  sin2  <px  +  b2  cos2  <pr  1 


(3.7) 


where  the  condition  that  |y>j|  <  tan  1(6/ y/d2  —  a2),  when  d  >  a,  and  x— tan  1(b/y/d2  —  a2)  < 


PT  <  « +tan-1(6/\/d2— a2),  when  d  <  —a,  is  imposed  to  insure  the  intersection  exists.  Sim¬ 
ilarly.  the  intersection  of  the  receiver  line-of-sight  with  (3.2)  can  be  obtained  by  substituting 
y  =  tan  pR(x  —  Lg)  into  (3.2);  thus,  we  have 

a2Lg  sin2  (pR  +  b2d  cos2  (pR  +  ab  cos  p>R\J[^  —  (d  —  Z,fl)2]  sin2  +  62  cos2 


a2  sin2  tpR  +  b2  cos2  <pR 


(3.8) 
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and 


y  = 


b2(d  -  La)  sm^RCOS  pR  +  ab  turn?  r  cos  (pRyj\a*-  (d  -  Lg)2]s\n2  ^r  +  62  cos2  kpR  /o  m 

— - - — •  ~  :  5  r~L2  i  ’ 

a2  sin  <pr  +  bl  cos2  <?r 


where  the  condition  that  |p*|  <  tan-1(6/^/(d  -  Tg)2  -  a2),  when  (d  -  Lg)  >  a,  and  tt  - 
tan"1  (&/>/(<* -£f)2-o2)  <  V*r  <  *  +  tan ^(ft/^d  -  Ta)2  -  a2),  when  (d  -  I3)  <  -a,  is 
imposed  to  insure  the  intersection  exists. 

The  derivations  in  this  section  are  based  on  the  assumption  of  flat  earth.  This  assumption 
is  applicable  in  practice  if  the  illuminated  region  is  small  such  that  it  can  be  approximated 
by  a  plane.  In  the  next  section  we  use  the  results  presented  here  to  analyze  the  ground 
clutter  Doppler  spread  in  airborne  bistatic  radar. 


3.2  Doppler  spread  in  ground  clutter 

The  previous  section  gives  the  formula  of  the  constant  range  contour  on  the  ground,  assuming 
flat  earth.  The  illuminated  portion  of  this  countour  defines  the  cross-range  dimension  of  the 
range  cell.  In  the  following  we  derive  the  Doppler  relationships  along  the  aforementioned 
contour.  It  will  be  seen  that  the  Doppler  of  the  ground  return  is  both  range  and  angle 
dependent.  The  ground  clutter  Doppler  spread  is  then  analyzed. 

Assume  that  the  velocity  of  the  transmitter  is  Vj,  in  the  direction  (cos  ~pvT  cos  0vT,  sin  <pvT  cos  9yT 
-  sin  0\-T ) .  and  the  velocity  of  the  receiver  is  Vr,  in  the  direction  (cos  <pvR  cos  9\-R,  sin  fvR  cos  9yR, 
-sin0vfl).  The  resulting  Doppler  at  a  point  on  the  ground,  (x,y,-hr),  can  be  written  as 

Id  =  --(*  cos  tpvT  cos  0vT  +  y  sin  <pvT  cos  9yT  +  hr  sin  0vT ) 

XRt 

+-^-[(x  —  Lg)  cos  ipyR  cos  9vr  +  y  sin  <f>vR  cos  0vR  +  Jir  sin  9vR),  (3.10) 

XRr 

where  X  is  the  wavelength.  Figure  3.1  to  Figure  3.3  show  the  constant  Doppler  contours  using 
the  same  scenario  as  given  in  [16]  but  with  the  restriction  of  near  zero  altitude  removed.  The 
values  used  are  hr  =  O.IjL,  Jir  =  0.2 L,  Vi  —  lOOm/s,  Vr  =  250m/s,  6vT  —  9yR  =  0,  and 
A  =  0.03m.  The  unit  of  range  in  the  figures  is  normalized  with  the  baseline  length.  Fig. 

3.1  depicts  the  case  when  ipvT  =  <PvR  =  tt/2,  while  Figure  3.2  depicts  the  case  ipvT  = 
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"/2*  r»>  =  ott/4.  It  can  be  seen  from  these  figures  that  the  Doppler  shift  varies  along 
the  ground  projection  of  either  the  transmitter  or  the  receiver  line-of-sight.  Consequently, 
collecting  secondary  data  along  either  direction  requires  different  Doppler  compensations 
for  each  secondary  range  cell  to  reconcile  with  the  clutter  Doppler  shift  at  the  test  cell. 
Nevertheless,  the  ground  clutter  Doppler  spread  may  still  depend  on  the  range  and  angle. 
In  the  following  we  examine  the  Doppler  spread  at  given  transmitter  and  receiver  azimuth 
angles. 

The  intersection  of  the  transmitter  line-of-sight  and  the  constant  range  contour  on  the 
ground  was  given  by  (3.6)  and  (3.7).  Substitution  of  these  expressions  into  (3.10)  gives  the 
Doppler  shift  at  that  point.  In  other  words,  Doppler  shift  at  the  points  on  (3.2)  is  expressed 
as  a  function  of  the  ranges  and  ipT-  Thus,  writing  the  Doppler  shift  as  fD(R,hT ,  hR,?T),  it 
is  easily  seen  that  the  maximum  change  in  Doppler  shift  over  the  constant  range  contour, 
within  the  transmitting  beam,  can  be  obtained  by 

A/Dr=  max  fD{R,  hT,hR,<pT)  -  min  fD(R,hT,hR,<pT),  (3.11) 

vt€DBt 

where  DRt  is  the  interval  of  ipx  of  the  transmitting  beam.  In  fact,  A  foT  is  the  ground  clutter 
Doppler  spread  when  the  cross-range  dimension  of  the  transmitting  beam  is  the  determining 
factor  of  the  cell  length.  Likewise,  substitution  of  (3.8)  and  (3.9)  into  (3.10)  may  be  used 
to  calculate  the  clutter  Doppler  spread  when  the  cell  length  is  determined  by  the  receiving 
beam: 

A /dr=  max  fD{R,hT,hR,<pR)  -  min  fD(R,hT,hR,pR),  (3.12) 

where  DRr  is  the  interval  of  <pR  of  the  receiving  beam. 

In  the  following  we  give  an  example  to  illustrate  how  the  Doppler  spread  varies  along 
a  given  direction  of  the  transmitting  beam;  similar  results  can  be  expected  along  a  given 
receiving  beam  direction.  The  beam  width  is  assumed  to  be  4 deg  and  the  geometry  in  Figs. 
3.1  and  3.2  are  used.  Figures  3.3  and  3.4  show  the  Doppler  spread  for  <px  =  Odeg  and 
yT  —  45 deg.  respectively. 

In  this  section  we  derived  the  Doppler  shift  along  a  constant  range  contour  on  the  ground 
for  an  airborne  bistatic  radar.  It  was  shown  that  the  Doppler  shift  is  both  range  and 
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Hertz 


Doppler  Spread 


Figure  3.3:  Clutter  Doppler  spread  along  the  direction  of  <pr  =  0  d 
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angle  dependent.  The  clutter  Doppler  spread  along  a  given  transmitting  or  receiving  beam 
direction  was  then  analyzed.  The  results  were  illustrated  with  an  example.  In  the  next 
chapter  we  examine  how  the  variation  in  the  clutter  spread  affects  the  estimation  of  the 
clutter  covariance  matrix. 
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Chapter  4 

A  Preliminary  Study  of  Space- Time 
Processing  for  Multistatic  Radar 
System 


In  this  chapter  we  present  some  results  from  a  preliminary  study  of  the  space-time  processing 
for  multistatic  system.  This  analysis  will  provide  us  with  a  basis  for  further  development. 
Attention  is  restricted  to  the  single  transmitter  case.  The  extension  to  a  more  generalized 
system  is  straightforward  and  involves  the  same  concept. 


Optimum  processing  and  the  steering  signal 

Suppose  the  multistatic  system  consists  of  one  transmitter  and  M  receivers.  The  returned 
signal  vector  is  now  given  as  s  =  [sfSj  where  Sm,m  =  1,2 is  the  signal 

sample  vector  obtained  at  the  mth  receiver.  The  vectors  of  the  noise  and  the  weights  of  the 
optimum  filter  can  be  formed  in  a  similar  manner. 

On  the  one  hand,  using  the  Schwarz  inequality,  the  same  conclusion  can  be  drawn  that  the 
optimum  filer  weights  are  of  the  form  given  by  (1.2);  on  the  other  hand,  a  direct  realization 
may  be  prohibitive  due  to  the  heavy  computational  requirement.  It  is  of  interest  to  find  out 
the  possibility  of  further  simplifying  (1.2). 

Assume  the  interferences  from  different  receivers  have  zero  mean  and  are  uncorrelated 
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such  that 


Rn  =  E hhH  = 


Rni  o 

0  Rn2 


R 


"A i  J 


(4.1) 


where  Rk.k  =  1.2,**-,  \f  are  the  noise  covariance  matrix  of  the  Ath  receiver  and  0  is  a  zero 
matrix.  It  follows  that 

k:  o  •  o  i 

«•-  ?  *  ... 

.  0  KL. 

Substituting  the  above  into  (1.2),  we  see  that  the  optimum  weights  of  each  receiver  can  be 
generated  independently. 

The  use  of  (4.2)  relies  on  whether  it  is  appropriate  to  assume  that  the  noise  at  one  receiver 
is  independent  from  that  from  the  others.  In  fact,  such  assumption  is  applicable  when  the 
thermal  noise  is  dominant,  but  it  may  not  be  suitable  when  clutter  is  the  prevailing  factor  of 
the  interference.  This  is  one  of  the  problem  to  be  solved  in  developing  space-time  processing 
for  multistatic  radar. 

There  is  another  problem  generally  overlooked,  namely,  the  employment  of  equally  spaced 
array  sensors.  This  configuration  may  not  cause  difference  in  the  monostatic  case,  where 
geometry  does  not  change  the  clutter  characteristics.  However,  the  analysis  in  Chapter 
2  indicates  that  the  bistatic  geometry  affects  the  observed  clutter  characteristics.  Conse¬ 
quently,  array  sensors  that  are  not  equally  spaced  may  be  more  suitable  for  the  application 
of  space-time  processing  in  the  bistatic/multistatic  case.  This  problem  is  to  be  analyzed. 
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Chapter  5 
Discussion 


In  this  report,  we  presented  an  analysis  of  the  clutter  characteristics  as  observed  from  a 
radar  platform  of  either  the  monostatic  or  the  bistatic  configuration,  airborne  or  stationary. 
This  study  shows  the  effects  of  system  parameters,  particularly  the  system  geometry,  on  the 
observed  clutter  characteristics.  However,  due  to  the  lack  of  a  suitable  model,  variation  of 
the  clutter  strength  in  the  out-of-plane  geometry  in  the  bistatic  radar  was  not  included.  This 
problem  will  be  considered  further  in  the  next  phase  of  this  research  effort. 

A  preliminary  study  of  adaptive  space-time  processing  for  multistatic  active  radar  can  be 
found  in  this  report  as  well.  This  analysis  show's  where  the  problems  lie  and  provides  us  with 
a  basis  for  further  development.  It  is  seen  that  there  are  two  major  factors  that  need  to  be 
investigated  further  in  the  development  of  space-time  processing  for  the  multistatic  system: 
(1)  the  correlation  of  the  observed  interference  among  the  receivers;  and  (2)  the  effects  of 
system  geometry  on  the  observed  reverberation  characteristics. 

As  the  analysis  shown,  the  problem  of  multistatic  space-time  processing  reduces  to  the 
monostatic  one  when  the  interferences  are  uncorrelated.  However,  this  may  not  be  the  case, 
e.g..  when  multipath  occurs.  Then,  it  is  of  interest  to  know  if  the  a  priori  knowledge  of  the 
system  geometry  can  be  used  to  reduce  the  computational  burden.  It  was  also  seen  from  the 
study  of  clutter  characteristics  that  the  observed  clutter  characteristics  does  depend  on  the 
system  geometry,  stationary.  Then,  the  question  is  how  it  affects  the  system  performance 
and  if  it  can  be  rectified.  Problem  also  arises  when  the  actual  signal  angle  of  arrival  does 
not  coincide  with  the  estimated  one,  or  when  the  actual  geometry  of  the  array  differs  from 
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the  one  assumed.  When  such  errors  occur,  the  actual  signal  may  be  tuned  out,  instead  of 
the  intended  interference. 

The  aforementioned  issues  are  essential  to  the  bistatic/multistatic  space-time  processing. 
They  need  to  be  addressed  before  an  effective  system  can  be  implemented.  The  solution  to 
them  is  the  goal  of  the  next  phase  of  this  research  effort. 
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Appendix  A 

Some  Results  on  Bistatic  Geometry 


Here  we  give  a  detailed  derivation  of  the  geometry  relationships  for  bistatic  radar.  Some 
special  cases  are  considered  for  illustration. 

A.l  Relationships  among  the  distance  parameters 

From  (2.43).  (2.54)  and  (2.55),  it  can  be  shown  that 

Rt  =  k-R’  +  B)  +  +  -  A1)}  (A.l) 

for  R  >  hr  +  \J L2  cos2  9rt  +  h2R,  and 

Rt  =  [R{R2  +  B)±  A^(R2  +  B)2  -  4h2T(R2  -  A2)}/  [l{R2  -  A2)]  (A.2) 

for  R  <  hr+  \J  L2  cos2  9rt  +  h2R.  Using  the  definition  that  R  =  Rt  +  Rr~  it  follows  that 
Rr  =  -  2 A1  -  B)  -  Aj(R?  +  Bf  -  thUR?  -  A1)]/ [2(tf  -  A2)]  (A.3) 

for  R  >  I>t  4-  \J L2  cos2  Ort  +  h\,  and 

Rr  =  [fl(fi2  -  2A2  -B)T  Ay/i^TW1 4/t^2-^2)]/  [2(^2  -  A2)]  (A.4) 

for  R  <  Iit  +  \J L2  cos2  9ftr  +  hR. 

Differentiating  (2.55),  we  have 

^  =  (-2  AR(A2  +  R)7(R2  +  R)2-  4/4(R2-A2) 
dR  l 

±[r?~  3 A2R*  -  (4 A2B  +  B2  -  4 A2h2T)R2  -  ( A2B 2  +  4 A*h2T)] } 

•  [2(i?2  -  A2)2sf[R2  +  B)2  -  4/4(R2  -  A2)]1 .  (A.5) 
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The  subtraction  in  (A. 5)  only  occurs  when  R  <  hr  +  L2  cos2  Brt  +  h2R.  Using  the  results 
for  Rt  and  Rr,  we  have 

=  ±2RtRr  [(tf  +  Bf  -  ih\(R?  -  X2)]  ',/2  ( A.6) 

the  negative  sign  only  occurs  when  R  <  fix  +  \Jl2  cos2  6rt  +  h2R. 

Special  cases 

Clearly,  the  parameters  for  the  monostatic  case  can  be  obtained  by  letting  L  =  0,  and  they 
are  not  derived  here.  When  hr  —  0, 


Rt  =  ( R 2  -  L2)/  [2(i?  -  L  cos  0RT  sin  <^r)] ,  (A.7) 

Rr  =  [R2  -  2 L2  cos2  0rt  sin2  ipr  +  hR  +  L2  cos2  6rt\/2(R  —  L  cos  Ort  sin  v?x),  (A.8) 

Rt9  =  Rt,  (A.9) 


A. 2  Width  of  the  range  ring  and  the  distance  to  the 
receiver 


From  the  definition  that  R  =  Rj  +  Rr,  we  have 

<1Rr _  dRj 

~dR  =  ~~dR‘ 

It  follows  that 

dRR  _  /  dRT\  dR 
Rr  ~  [  dR  )  Rr' 

Using  (2.54),  (2.67)  and  (A.3)  in  (A.  12),  we  have 


(A.ll) 

(A.12) 


<IRr  _  J _ 2(1  -  A2/fi2) _ 

Rr  ~  \  1  -  2 A2/R2  -  B/R 2  -  (A/R)y/{1  +  bJr2)2  -  4A£/i?2(l  -  A2/#2) 

(A/m  +  B/&)  +  ^(1  +  W&y  -  4fej/fi2(l  -  AV~^))  dR 
(1  -  A2/R2)J(  1  +  fi/#2)2  -  4/4/fl2(l  -  A2/#2)  J 
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Invoking  the  assumptions  that  the  square  of  the  ratios  hT/R*  hR/R  and  L/R  are  much 
smaller  than  unity,  the  above  relationship  reduces  to 


dRR  1  dR 
Rr  “  1  -A/R  R' 


(A.14) 


Also,  following  the  assumptions,  the  second  and  higher  powers  of  \A/R\  are  much  smaller 


than  unity;  thus. 


Since  \A/R\  is  smaller  than  unity,  we  conclude  that 


dR 
R  ‘ 


(A.15) 


ARr/Rr  <  1  when  A R/R  <  1  within 


one  range  ring. 
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